NONLINEAR ACCELERATOR PROBLEMS VIA WAVELETS: 
8. INVARIANT BASES, LOOPS AND KAM 
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In this series of eight papers we present the applications of 
methods from wavelet analysis to polynomial approxima- 
tions for a number of accelerator physics problems. In this 
part we consider variational wavelet approach for loops, in- 
variant bases on semidirect product, KAM calculation via 
FWT. 



1 INTRODUCTION 

This is the eighth part of our eight presentations in which 
we consider applications of methods from wavelet analy- 
sis to nonlinear accelerator physics problems. This is a 
continuation of our results from [l]-[8], in which we con- 
sidered the applications of a number of analytical methods 
from nonlinear (local) Fourier analysis, or wavelet analy- 
sis, to nonlinear accelerator physics problems both general 
and with additional structures (Hamiltonian, symplectic or 
quasicomplex), chaotic, quasiclassical, quantum. Wavelet 
analysis is a relatively novel set of mathematical methods, 
which gives us a possibility to work with well-localized 
bases in functional spaces and with the general type of 
operators (differential, integral, pseudodifferential) in such 
bases. In contrast with parts 1-4 in parts 5-8 we try to 
take into account before using power analytical approaches 
underlying algebraical, geometrical, topological structures 
related to kinematical, dynamical and hidden symmetry of 
physical problems. In section 2 we consider wavelet ap- 
proach for calculation of Arnold- Weinstein curves (closed 
loops) in Floer variational approach. In section 3 we con- 
sider the applications of orbit technique for constructing 
different types of invariant wavelet bases in the particular 
case of affine Galilei group with the semiproduct structure. 
In section 4 we consider applications of very useful fast 
wavelet transform (FWT) technique (part 6) to calculations 
in KAM theory (symplectic scale of spaces). This method 
gives maximally sparse representation of (differential) op- 
erator that allows us to take into account contribution from 
each level of resolution. 



2 FLOER APPROACH FOR CLOSED 
LOOPS 

Now we consider the generalization of wavelet variational 
approach to the symplectic invariant calculation of closed 
loops in Hamiltonian systems [9]. As we demonstrated 
in [3]-[4] we have the parametrization of our solution by 
some reduced algebraical problem but in contrast to the 
cases from parts 1-4, where the solution is parametrized by 
construction based on scalar refinement equation, in sym- 
plectic case we have parametrization of the solution by ma- 
trix problems - Quadratic Mirror Filters equations. Now 
we consider a different approach. Let(M, uS) be a compact 
symplectic manifold of dimension 2n, u> is a closed 2-form 
(nondegenerate) on M which induces an isomorphism 
T*M — ► TM. Thus every smooth time-dependent Hamil- 
tonian H : R x M — » R corresponds to a time-dependent 
Hamiltonian vector field X H ■ R x M — * TM defined 
by u(X H (t,x),£) = -d x H(t,x)£ for £ e T X M. Let H 
(and Xh ) is periodic in time: H (t + T, x) = H(t, x) and 
consider corresponding Hamiltonian differential equation 
on M: x(t) = Xjj(t,x(t)) The solutions x(t) determine 
a 1-parameter family of diffeomorphisms ipt S Diff(M) 
satisfying i/j t (x(0)) = x(t). These diffeomorphisms are 
symplectic: uj = ipfUi. Let L = LtM be the space of 
contractible loops in M which are represented by smooth 
curves 7 : R — > M satisfying j(t + T) = j(t). Then 
the contractible T-periodic solutions can be characterized 
as the critical points of the functional S = St ■ L — > R: 
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where D C C be a closed unit disc and u : D — > M is 
a smooth function, which on boundary agrees with 7, i.e. 
u(exp{27ri9}) = j(QT). Because [ui], the cohomology 
class of uj, vanishes then St (7) is independent of choice 
of u. Tangent space T 7 L is the space of vector fields £ S 
C°°(7*TA/) along 7 satisfying f (i + T) = £(t). Then we 
have for the 1-form df : TL — > R 



dSr( 7 )£ 
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and the critical points of S are contractible loops in L 
which satisfy the Hamiltonian equations. Thus the criti- 
cal points are precisely the required T-periodic solutions. 



To describe the gradient of S we choose a on almost com- 
plex structure on M which is compatible with oj. This is an 
endomorphism J £ C°°(End(TM)) satisfying J 2 = -I 
such that g(£,rj) = J(x)rf), £,77 € T X M defines a 
Riemannian metric on M. The Hamiltonian vector field is 
then represented by Xn(t, x) = J(x)W H(t, x), where V 
denotes the gradient w.r.t. the x-variable using the metric. 
Moreover the gradient of S w.r.t. the induced metric on 
L is given by grad5(7) = .7(7)7 + Vff(t,7), 7 £ L. 
Studying the critical points of S is confronted with the 
well-known difficulty that the variational integral is nei- 
ther bounded from below nor from above. Moreover, at 
every possible critical point the Hessian of / has an infinite 
dimensional positive and an infinite dimensional negative 
subspaces, so the standard Morse theory is not applicable. 
The additional problem is that the gradient vector field on 
the loop space L: d7/ds = — grad/(7) does not define 
a well posed Cauchy problem. But Floer [9] found a way 
to analyse the space j\A of bounded solutions consisting of 
the critical points together with their connecting orbits. He 
used a combination of variational approach and Gromov's 
elliptic technique. A gradient flow line of / is a smooth 
solution u : R — > M of the partial differential equation 

nit nn 

^ + J(u)^ + VH(t,u) = 0, (3) 
os at 

which satisfies u(s, t + T) = u(s, t). The key point is to 
consider ^ not as the flow on the loop space but as an el- 
liptic boundary value problem. It should be noted that (^J) 
is a generalization of equation for Gromov's pseudoholo- 
morphic curves (correspond to the case WH = in 
Let M.t = A4t(H, J) the space of bounded solutions of 
([3]), i.e. the space of smooth functions u : C/iTZ — > M, 
which are contractible, solve equation (|j) and have finite 
energy flow: 

•"*•>- 5/ T(l^ 2 + !§F- X '*»>l 2 ) dtd »- 

(4) 

For every u £ Mt there exists a pair x, y of contractible 
T-periodic solutions, such that u is a connecting orbit from 
y to x: linis^oo u(s,t) = y(t), lim s ^ +00 = x(t). Then 
our approach from preceding parts, which we may apply or 
on the level of standard boundary problem or on the level of 
variational approach and representation of operators (in our 
case, J and V) according to part 6(FWT technique) lead us 
to wavelet representation of closed loops. 

3 CONTINUOUS WAVELET 
TRANSFORM. BASES FOR 
SOLUTIONS. 

When we take into account the Hamiltonian or Lagrangian 
structures from part 7 we need to consider generalized 
wavelets, which allow us to consider the corresponding 
structures instead of compactly supported wavelet repre- 
sentation from parts 1-4. We consider an important par- 
ticular case of constructions from part 7: affine relativity 



group (relativity group combined with dilations) — affine 
Galilei group in n-dimensions. So, we have combination of 
Galilei group with independent space and time dilations: 
Gaff = G rn tx D 2 , where D 2 = (R+) 2 ~ R 2 , G m 
is extended Galilei group corresponding to mass param- 
eter m > (G a ff is noncentral extension of G ex D 2 
by R, where G is usual Galilei group). Generic element 
of G a ff is g = ($, b n , b; v ; R, ao, a), where $ e R is 
the extension parameter in G m , bo £ R, b £ R™ are 
the time and space translations, v £ R" is the boost pa- 
rameter, R £ SO(n) is a rotation and ao,a £ R+ are 
time and space dilations. The actions of g on space-time 
is then x 1— > aRx + ciovt + b, t 1— > a^t + bo, where 
x = (xi,x 2 , ...,x n ). It should be noted that D 2 acts non- 
trivially on G m . Space-time wavelets associated to G a ff 
corresponds to unitary irreducible representation of spin 
zero. It may be obtained via orbit method. The Hilbert 
space is TC — i 2 (R" x R, dkduj), k — (fc l5 k n ), where 
R™ x R may be identified with usual Minkowski space and 
we have for representation: 

(U(g)^)(k, w) = y/a a n expi(m§ + kb- uib )^(k', oj'), 

(5) 

with k' = aR^ 1 (k + mv), oj' = ao(u> — kv — ^mv 2 ), 
m' = (a 2 /ao)m. Mass m is a coordinate in the dual of 
the Lie algebra and these relations are a part of coadjoint 
action of G a ff- This representation is unitary and irre- 
ducible but not square integrable. So, we need to consider 
reduction to the corresponding quotients X = G/H. We 
consider the case in which H={phase changes $ and space 
dilations a}. Then the space X = G/H is parametrized 
by points x — (bo, b; v; R; ao). There is a dense set of vec- 
tors rj £ H admissible mod(H, 073), where 073 is the corre- 
sponding section. We have a two-parameter family of func- 
tions /^(dilations): f3(x) = (fio + Ajao) 1 / 2 , Ao,^o € R. 
Then any admissible vector 77 generates a tight frame of 
Galilean wavelets 

»fe(s)(*,w) = v/«o(mo + Aoao^/V^-^Vfc', w'), 

(6) 

with k' — (/i + \oa) 1 ^ 2 Rr 1 {k + mv), oj' — a (uj — 
kv — mv 2 /2). The simplest examples of admissible vec- 
tors (corresponding to usual Galilei case) are Gaussian 
vector: rj(k) ~ exp(— k 2 /2mu) and binomial vector: 
7]{k) ~ (1 + k 2 j2mu)^ a l 2 , a > 1/2, where u is a 
kind of internal energy. When we impose the relation 
ao = a 2 then we have the restriction to the Galilei- 
Schrodinger group G s — G m ex D s , where D s is the 
one-dimensional subgroup of D 2 . G s is a natural in- 
variance group of both the Schrodinger equation and the 
heat equation. The restriction to G s of the represen- 
tation (29) splits into the direct sum of two irreducible 
ones U = U+ U- corresponding to the decomposi- 
tion i 2 (R" x R,dfcdw) = H+ 8 H-, where H± = 
L 2 (D±,dkduj = {ip £ L 2 (R n x R,dfcdw), ip(k,w) = 
for oj + k 2 /2m = 0}. These two subspaces are the ana- 
logues of usual Hardy spaces on R, i.e. the subspaces of 
(anti)progressive wavelets (see also below, part III A). The 



two representation U± are square integrable modulo the 
center. There is a dense set of admissible vectors rj, and 
each of them generates a set of CS of Gilmore-Perelomov 
type. Typical wavelets of this kind are: the Schrodinger- 
Marr wavelet: r](x,t) = (id t + A/2m)e~ (;r2+i2) / 2 , the 
Schrodinger-Cauchy wavelet: ip(x,t) = (idt + A/2m) x 
(t + i) YYj—i (xj + i) ■ So, in the same way we can con- 
struct different invariant bases with explicit manifestation 
of underlying symmetry for solving Hamiltonian or La- 
grangian equations. 

4 SYMPLECTIC HILBERT SCALES VIA 
WAVELETS 

We can solve many important dynamical problems such 
that KAM perturbations, spread of energy to higher modes, 
weak turbulence, growths of solutions of Hamiltonian 
equations only if we consider scales of spaces instead of 
one functional space. For Hamiltonian system and their 
perturbations for which we need take into account under- 
lying symplectic structure we need to consider symplectic 
scales of spaces. So, if ii(t) = JW K (u(t)) is Hamilto- 
nian equation we need wavelet description of symplectic or 
quasicomplex structure on the level of functional spaces. It 
is very important that according to [12] Hilbert basis is in 
the same time a Darboux basis to corresponding symplec- 
tic structure. We need to provide Hilbert scale {Z s } with 
symplectic structure [12]. All what we need is the follow- 
ing. J is a linear operator, J : Z M — * Z^, J(Z oa ) = Z m , 
where Z M = DZ S . J determines an isomorphism of scale 
{Z s } of order dj > 0. The operator J with domain of def- 
inition Zoo is antisymmetric in Z: < Jzi,Z2 >z= — < 
zi,Jz2 >z,zi,Z2 G Zoo. Then the triple {Z,{Z s \s 6 
R}, a =< Jdz,dz >} is symplectic Hilbert scale. So, 
we may consider any dynamical Hamiltonian problem on 
functional level. As an example, for KdV equation we 
have Z s = {u(x) € H s (T 1 )^ u(x)dx = 0}, s G R, 
J = d/dx, is isomorphism of the scale of order one, 
J = — (J) -1 is isomorphism of order —1. According to 
[13] general functional spaces and scales of spaces such 
as Holder-Zygmund, Triebel-Lizorkin and Sobolev can be 
characterized through wavelet coefficients or wavelet trans- 
forms. As a rule, the faster the wavelet coefficients decay, 
the more the analyzed function is regular [13]. Most im- 
portant for us example is the scale of Sobolev spaces. Let 
Hk{R n ) is the Hilbert space of all distributions with finite 
norm \\s\\ 2 Hk{Rn) = /dtfl + \H\ 2 ) k / 2 \m\ 2 ■ Let us con- 
sider wavelet transform 



W g f(b,a) 



dx — q 

Rn a 



b 6 R n , a > 0, w.r.t. analyzing wavelet g, which is 
strictly admissible, i.e. C 9 , s = J^°(da/a)\g(ak)\ 2 < 00. 
Then there is a c > 1 such that 



W g s(b,a)\ 2 <c\\sf Hk{Rny 

This shows that localization of the wavelet coefficients at 
small scale is linked to local regularity. So, we need repre- 
sentation for differential operator (J in our case) in wavelet 
basis. We consider it by means of the methods from part 6. 

We are very grateful to M. Cornacchia (SLAC), W. Her- 
rmannsfeldt (SLAC) Mrs. J. Kono (LBL) and M. Laraneta 
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